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We demonstrate the existence of metallic spin density waves (SDWs) in the Kondo lattice model 
on a square lattice for a wide range of parameters by means of real space dynamical mean field theory. 
In these SDWs, the spin polarization as well as the charge density depend on the lattice site and are 
modulated along one direction of the square lattice. We show that within this phase of metallic SDWs 
the Fermi surface changes from small to large, when the coupling strength is increased. Furthermore, 
the transition between the large Fermi-surface SDW phase and the paramagnetic phase is of second 
order, while the transition between the small Fermi-surface SDW phase and the paramagnetic phase 
is of first order. A local quantum critical point is thus avoided in our calculations by undergoing a 
first order phase transition. 


PACS numbers: 71.10.Fd, 71.27.+a, 75.30.Fv 

I. INTRODUCTION 

Strongly correlated materials are of special interest in 
condensed matter physics because they often exhibit re¬ 
markable phenomena which cannot be explained by sin¬ 
gle particle theory. A particular class of strongly cor¬ 
related materials are heavy fermion systems, in which 
quantum criticality accompanied by non-Fermi liquid be¬ 
havior or unconventional superconductivity can often be 
observed.Heavy fermion materials include strongly 
interacting /-electrons which are hybridized with conduc¬ 
tion 5 pd-bands. The main driving forces in these materi¬ 
als are the Kondo screening, which favors a paramagnetic 
state, and the Ruderman-Kasuya-Kittel-Yosida (RKKY) 
interaction, which favors a magnetically long-range or¬ 
dered state. Especially, when both effects are of com¬ 
parable strength, fascinating many-particle phenomena 
such as quantum criticality and non-Fermi-liquid physics 
can arise. 

It has been found that quantum phase transitions 
in heavy fermion systems can be categorized into two 
classes: (i) Quantum fluctuations at the phase tran¬ 
sition arise solely due to the magnetic order and the 
Fermi surface behaves smoothly across the transition, 
which corresponds to the Hertz-Millis-Moriya type of 
quantum critical point (QCP),^^^®^ and (ii) exactly at 
the quantum phase transition between magnetism and 
paramagnetism, not only the magnetic order vanishes, 
but also the Kondo screening vanishes and the Fermi 
surface changes ab rupt ly, which is referred to as local 
quantum criticality.^^^ While CeNi2Ge2 is an example 
for the Hertz-Millis-Moriya type of QCP,^ there is evi¬ 
dence that the QCP observed in YbRh2Si2 corresponds 
to a local quantum critical point.^^^^^ These experimen¬ 
tal and theoretical results on quantum phase transi¬ 
tions in heavy fermion systems have been summarized 
in the global phase diagranP^^KSElHllI which includes 
small/large Fermi-surface antiferromagnetic phases, and 
small/large Fermi-surface paramagnetic phases. 

A commonly used theoretical model to describe heavy 


fermion systems is the Kondo lattice model, in which 
effective local moments originating from the interact¬ 
ing /-electrons are coupled to conduction electrons. The 
Kondo lattice model read^^^^^ 

H = t Y. cIc^^ + jY^SA (1) 

<i,j>(T i 

~ ^icFmP^rnCrn^icFn^ 

where c[^ creates an electron with spin-direction cr on lat¬ 
tice site i. The first term in Eq.Q describes the hopping 
of conduction electrons, and the second term describes 
the coupling of the conduction electrons to the localized 
spins with interaction strength J. Throughout this pa¬ 
per we assume an antiferromagnetic coupling, J > 0 , 
and take the hopping constant t as unit of energy. All 
calculations are performed on a two-dimensional square 
lattice. 

In this paper we examine magnetism and quantum 
phase transitions in the Kondo lattice model by means 
of the dynamical mean field theory (DMFT) ,I^^M^ which 
has been extensively used for this purposeAlthough 
these previous studies include calculations for the an¬ 
tiferromagnetic state and the magnetic quantum phase 
transition, they were performed for the Neel state with 
commensurate ordering vector. Exactly at half filling for 
a particle-hole symmetric lattice, the Neel state is insu¬ 
lating. Thus, the behavior of the Eermi surface across 
the quantum critical point cannot be studied. Some of 
these previous studies also used the Neel state away from 
half filling within the metallic antiferromagnetic phase. 
However, away from half filling, one can expect that the 
Neel state becomes energetically unstable towards incom¬ 
mensurate magnetic states, as has also been found in 
many experimental systems by neutron scattering, e.g. 
for CeRu2Si2!^ CeCu2Si2}^ CeRhlns,!^ etc. A theoret¬ 
ical study using classical spins already proved the exis¬ 
tence of incommensurate SDWs.l^ However, due to the 
approximation of classical spins, many-particle effects, 
particularly the Kondo effect, are absent in Ref. 61 . 
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In this paper, we use the DMFT, and explicitly in¬ 
clude the possibility of incommensurate magnetic states 
in the Kondo lattice model by using large-scale real-space 
DMFT (RDMFT). We will focus on the analysis of metal¬ 
lic SDW states away from half filling, which have not 
been carefully studied yet, and analyze static and dynam¬ 
ical properties with particular emphasis on the behavior 
of the Fermi surface across the quantum phase transition. 

This paper is organized as follows: In the next sec¬ 
tion (II), we will shortly explain the methods used in 
this paper. This is followed by the section (III), in which 
we will summarize the calculated phase diagram and ex¬ 
plain the different types of SDWs observed in our calcu¬ 
lations. Static (IV), dynamical properties (V) and the 
phase transitions between different phases (VI) are ex¬ 
plained in more detail in the subsequent sections, before 
concluding the paper. 


II. METHOD 

We use large-scale real-space dynamical mean field the¬ 
ory (RDMFT) to study the Kondo lattice model, Eq. 0, 
defined on a two-dimensional square lattice. In RDMFT 
each lattice site of a finite cluster is mapped onto its own 
impurity model. This mapping is done by calculating the 
local Green’s function of each lattice site via 



Figure 1: Phase diagram of the Kondo lattice model as cal¬ 
culated by RDMFT. The SDW phase away from half filling 
is separated into a large Fermi-surface (FS) phase at strong 
coupling and a small Fermi-surface phase at weak coupling. 
A detailed explanation is given in the main text. 


“mixing” and frequently used in DMFT calculations. Re¬ 
cently, we have used similar calculations to study spin- 
and charge density waves in the Hubbard model.l^ 


G 


loc 



dhxdky • H 



( 2 ) 


III. PHASE DIAGRAM 


where is the hopping Hamiltonian of the chosen 

cluster and is the self-energy matrix. The momentum- 
dependence of arises thereby through the bound¬ 

ary conditions. For this study, we use (20x20)-, (60x 10)- 
and (30 x 30)-cluster with periodic boundary conditions. 
The self-energy is diagonal and momentum-independent 
within this approximation; correlation effects between 
different lattice sites are not included. However, the self¬ 
energy depends on the lattice site of the cluster, which 
makes it possible to describe SDW states. From the local 
Green’s function, a lattice-site-dependent hybridization 
is calculated via 


Aii{z) = z- 


Gloc,ii (^) 


(^) •) 


( 3 ) 


which determines the input for a single impurity Kondo 
model. These equations are reduced to the usual single¬ 
site DMFT, if the cluster consists of a single site. To 
obtain the self-energy of the Kondo impurity m odel, 
we use the numerical renormalization group (NRG),l^^*^ 
which can ca lculat e reliable self-energies at zero and finite 
temperature l^^l^^ l We then iterate Eqs. ([^ and ^ until 
the self-energy matrix becomes self-consistent. In order 
to improve the stability of the iterative calculations, we 
do not directly use the self-energy as calculated by the 
NRG, but we add 50% of the self-energy of the last it¬ 
eration, = 0.5E^ + 0.5E^_i, which is usually called 


Eigure summarizes our obtained results. Close to 
half filling, the physics of the Kondo lattice model is 
dominated by the interplay of the Kondo effect and the 
RKKY interaction, as described by the Doniach phase 
diagram.l^ At weak coupling, where the RKKY interac¬ 
tion dominates, we observe a magnetically ordered phase. 
At strong coupling, where the Kondo effect becomes dom¬ 
inating, a paramagnetic state is stabilized. 

Exactly at half filling the model exhibits perfect nest¬ 
ing and we observe the well-known transition from an 
antiferromagnetic Neel state at weak coupling to a para¬ 
magnetic Kondo insulator at strong coupling. This con¬ 
tinuous transition occurs approximately at Jc/t ^ 2.2 
within DMET for the square lattice. Notice that the 
Kondo lattice model on a square lattice remains insulat¬ 
ing at half filling for all interaction strengths, J > 0, al¬ 
though the origin of the insulating gap gradually changes. 
Our obtained critical interaction strength for the transi¬ 
tion between the Neel state and the Kondo insulator rea¬ 
sonably agrees with previous DMET calculations.^^ We 
can thus reproduce previous results with our RDMET 
calculations. 

Away from half filling, the homogeneous Neel state be¬ 
comes unstable and changes into a phase of SDWs, in 
which the polarization of the conduction electrons as well 
as the localized spins are lattice-site dependent. Eor cou¬ 
pling strengths larger than the hopping amplitude, this 
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Figure 2: Different magnetic states close to half filling: a) an¬ 
tiferromagnetic Neel state exactly at half filling {J — 2t). b) 
doped SDW/CDW at weak interaction strength (J = 0.8t). 
c) doped SDW/CDW close to the quantum critical phase 
transition (J = 1.8t). For each state a)-c), we show a false- 
color-plot of the charge density (top left), a false-color-plot 
of the electron polarization (top right), an intersection of the 
charge-density plot for a given y-siie (bottom left), and an 
intersection of the electron-polarization plot for a given y-site 
(bottom right). 


SDW phase becomes unstable upon doping towards the 
large Fermi-surface paramagnetic state. On the other 
hand, for weak coupling, J < t, SDWs can be observed 
for any conduction band filling. 

Even within this SDW phase, the modulation of the 
electron polarization and electron density depends on the 
model parameters. In Fig. we compile information 
about typical SDWs, which can be observed in the vicin¬ 
ity of half filling. For each type of SDW, a)-c), we show 
in the upper part of each panel false-color plots of the 
electron density (left plot) and the electron polarization 
(right plot) for a cluster of lattice sites. The lower plots 
of each panel show intersections of the above false-color 
plots. 

Panel a) shows the antiferromagnetic Neel state ex¬ 
actly at half filling. In the Neel state, the polarization of 
the electrons changes its sign between nearest neighbors. 
The electron density, on the other hand, is unity for all 
lattice sites. Increasing the interaction strength beyond 
J/t = 2.2, the antiferromagnetic state changes into the 
Kondo insulator (not shown in Fig. [^, where all lattices 
sites are occupied in average with a single electron and 
the polarization vanishes. 

Away from half filling, the perfect nesting, which ex¬ 
ists exactly at half filling, is lost, and the system can gain 
energy by modulating the electron polarization over the 
lattice sites. We have tried in our calculations different 
types of modulations, such as vertical SDWs and diag¬ 
onal SDWs. It turned out that the most stable type of 
SDW on the square lattice is the vertical SDW, as shown 
in panel b) and c) of Fig. The amplitude of the elec¬ 
tron polarization is modulated along one direction of the 
square lattice. Along the other direction, the polariza¬ 
tion shows a Neel-state-like ’AB’ oscillation. However, 
not only the polarization is modulated, but also the elec¬ 
tron density depends on the lattice site; there are stripes 
of high electron density and stripes of low electron den¬ 
sity. CDWs are well known in strongly correlated ma¬ 
terials, e.g. High-Tc Cuprates, and have been detected 
there by X-ray scattering.!^ Within the regions of high 
electron density, we generally observe a Neel-state-like 
behavior; the polarization changes its sign between near¬ 
est neighbors. The behavior in the low density regions, 
however, depends on the interaction parameter and the 
average occupation of the conduction electrons. For large 
interaction strengths close to the quantum critical transi¬ 
tion, we observe that the polarization vanishes in the low 
electron-density regions, see Fig. panel c). The Kondo 
effect becomes strong in these regions so that the local¬ 
ized moments are screened by the Kondo effect instead 
of forming a magnetically ordered state. On the other 
hand, for weak interaction strength, J < t, and espe¬ 
cially for an average electron filling n < 0.9, we observe 
that the low electron-density regions are still magneti¬ 
cally polarized. We find that nearest neighbor sites form 
ferromagnetic bonds along one axis, see Fig. panel b). 
For weak interaction strengths, the localized moments 
cannot be screened by the Kondo effect, which is expo- 
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Figure 3: Upper panel: Electron polarization in the high den¬ 
sity regions (corresponding to the maximal polarized elec¬ 
trons) for different coupling strengths and conduction band 
Filings. The red lines correspond to contour lines at \n^ — 
ni\ = (0; 0.1; 0.2). Lower panel: The difference of the elec¬ 
tron density between high and low density regions. The red 
lines correspond to contour lines at jnhigh — '/T'IowI = (0; 0.03; 
0.06; 0.09; 0.12). 


nentially weak. Here the magnetic state is completely 
determined by the RKKY interaction. The localized mo¬ 
ments are nearly fully polarized and form also inside the 
low electron-density region a magnetic state. Further¬ 
more, the modulation of the polarization cannot be de¬ 
scribed anymore by a simple sine-wave, but a more com¬ 
plex polarization pattern is realized. Due to a changed 
Fermi vector in this region, the system can gain energy by 
combining ferromagnetic bonds with antiferromagnetic 
ones. Experimentally, similar combinations of antiferro¬ 
magnetic and ferro magn etic correlations have been ob¬ 
served in YbRh 2 Si 2 ! ^^ * ^ * 


Eigure 4: Strongest Eourier mode of the charge density (up¬ 
per panel) and spin polarization (lower panel) for different 
conduction band Filings. The interaction value is Jjt — 0.8, 
for which the SDW phase exists until (n) 0.6. 


IV. STATIC PROPERTIES OF THE SDW 
PHASE 


An important point which was neglected in previous 
studies is that both the spin polarization and the elec¬ 
tron density are modulated away from half Filing. In 
Fig. 1^ we show the amplitude of the SDW (upper panel) 
and the amplitude of the CDW (lower panel) for differ¬ 
ent coupling strengths and conduction band Filings. We 
observe that the electron polarization within the SDW 
is largest for approximately J ^ t exactly at half Fil¬ 
ing. The spin polarization becomes weaker when dop¬ 
ing the system away from half Filing. The amplitude 
of the CDW, however, is strongest for J « t and con¬ 
duction band Filings (n) 0.9. Exactly at half Filing, 

the electron density is homogeneous. Furthermore, the 
CDW is coupled to the SDW; within the paramagnetic 
state, no charge-density modulation exists. The Kondo 
effect alone is not sufficient to stabilize a CDW state for 
these conduction band Filings. With increasing coupling 
strength, the CDW vanishes faster than the SDW. How¬ 
ever, this might be related to the fact that the SDW state 
can only be found very close to half Filing for strong in¬ 
teraction. In this parameter region the ground state be¬ 
comes more and more homogeneous, and the amplitude 
of the CDW is very small. 

The observed modulations of the spin polarization and 
the electron density are not independent of model param¬ 
eters. In order to obtain information about the wave¬ 
length of the modulations, we have performed a Fourier 
transformation of the obtained DMFT solutions. For this 
purpose, we have used (60 x 10)-clusters for our calcu¬ 
lations. The Fourier transformation for the occupation 
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and the polarization read 
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where we have subtracted the average filling of the lattice, 
naver = jj Ea;,y (%(*,?/)+%(*,?/))> ^nd N Corresponds to 
the number of lattice sites. The absolute values of these 
Fourier modes show distinct peaks at certain momenta in 
the Brillouin zone. We show these momenta in Fig. |^for 
J/t = 0.8. For this interaction strength, the SDW phase 
exists until close to quarter filling, where it changes into 
another magnetic phase via a first order transition. The 
wavelength does not depend on the coupling strength J 
within the weak coupling region of the phase diagram. 
If we neglect the parameter region of small hole doping, 
where the wavelength of the SDW exceeds the cluster 
size, the largest Fourier component in the SDW and the 
CDW can be fitted by a linear function in the number of 
holes. Exactly at half filling, we find a Neel state with 
modulation k = (tt, tt) without a CDW. Away from half 
filling this Neel state changes into an SDW with mod¬ 
ulation k = 7r(l — (h), 1) accompanied by a CDW with 


modulation k = 7r(2(h),0), where (h) = 1 — {n) corre¬ 
sponds to the number of holes in the lattice. For stronger 
interaction values, the SDW phase vanishes upon doping 
into the paramagnetic heavy fermion state for small hole 
doping. The wavelength in this parameter region is too 
long, as that we could accurately analyze the doping de¬ 
pendence of it. 


V. SPECTRAL FUNCTION AND FERMI 
SURFACE 

As shown in Fig. we have divided the SDW phase 
into two regions; one with large and one with small Fermi 
surface. We next want to show momentum-resolved 
spectral functions and the corresponding Fermi surfaces 
which justify this distinction. 

Figure shows spectral functions for the paramagnetic 
phase (left), the strong-coupling SDW phase (middle), 
and the weak-coupling SDW phase (right). The lower 
panels correspond to magnifications around the Fermi 
energy of the upper panels. The green dashed line marks 
the Fermi energy. In the weak-coupling SDW phase 
(right panel in Fig. |^, we only find slight modifications 
from the non-interacting energy momentum dispersion. 
These modifications mainly occur close to the Fermi en¬ 
ergy around (7r,0), (0,7r), and (7r/2,7r/2), where we ob¬ 
serve a suppression of the spectral weight. Compared 
to the non-interacting spectrum, additional weak bands 
occur around (0, 0) and (tt, tt) below the Fermi energy. 
Especially the band which is close to the Eermi energy 
at (tt, tt) is a remnant of the heavy fermion band, which 


however lies below the Eermi energy. The main contribu¬ 
tion to the spectral weight at the Eermi energy originates 
from the non-interacting dispersion at (7r,0), (0,7r), and 
(7r/2,7r/2). 

The spectral functions of the paramagnetic metallic 
phase (left panel in Eig. and the strong-coupling SDW 
phase (middle panel), differ substantially from the weak- 
coupling case. The bands, which cross the Eermi en¬ 
ergy in the non- and weakly-interacting case at (tt, 0) 
and (0,7r), are absent at the Eermi energy in this phase. 
While in the weakly-interacting SDW phase there is a 
band which connects the low energy parts {ujjt = —4) 
with the high energy parts {uj/t = 4), in the param¬ 
agnetic and the strong-coupling SDW state there is a 
gap above the Eermi energy separating the low and the 
high energy parts of the spectrum. The non-interacting 
bands, which approach the Eermi energy at (tt, 0), (0, tt), 
and (7r/2,7r/2) are bent towards each other and cross the 
Eermi energy close to (7r,7r). The spectrum of the para¬ 
magnetic state and the strong-coupling SDW state look 
thereby very similar to each other. The gap above the 
Eermi energy seems to be more pronounced in the SDW 
state. However, this is parameter dependent. Another 
general feature in the dispersion of the SDW state is that 
there is a band at (0, 0) close to the Eermi energy. This 
is a shadow band originating from the ordered state.l^^^^ 

Erom the spectral functions shown in Eig. we can see 
that the Eermi surface changes within the SDW phase. 
Eor a better understanding, we have compiled the cal¬ 
culated Eermi surfaces for a few parameters in Eig. 
Different rows show different transitions, which can be 
observed in the Kondo lattice model. The upper pan¬ 
els show the transition between the weak-coupling SDW 
phase and the paramagnetic phase at (n) ~ 0.8 for J/t = 
(0.6; 0.8; 1.0; 1.2; 1.4) (left-right). One can clearly observe 
a change in the Eermi surface topology at J/t = 1. While 
in the left two panels the Eermi surface is small including 
only conduction electrons, the Eermi surface in the right 
two panels is large, which includes conduction electrons 
and localized moments. A similar change in the Eermi 
surface topology can be observed when increasing the in¬ 
teraction strength within the SDW phase (middle panels 
of Eig. [^. The Eermi surface clearly changes in the sec¬ 
ond panel from right for J/t = 1.2 and (n) 0.9. Einally, 

the lower panels show the transition between the strong- 
coupling SDW and the paramagnetic phase for different 
conduction band fillings and J/t = 1.8. The magnetic 
order vanishes in the middle panel for (n) 0.93. How¬ 

ever, there is no abrupt change in the Eermi surface for 
strong coupling when the magnetic order vanishes. Sum¬ 
marizing these results, we find that the SDW phase away 
from half filling is divided into a phase with small Eermi 
surface at weak coupling and a phase with large Eermi 
surface at strong coupling. We thus find a Lifshitz tran¬ 
sition within the SDW phase, where the Eermi surface 
topology changes. 

Although the Eermi surface vanishes exactly at half 
filling due to the insulating nature of the system, a 
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Figure 5: Spectral functions for the paramagnetic state (left), the strong-coupling SDW (middle), and the weak-coupling SDW 
state (right). The lower panel is a magnihcation of the upper panel around the Fermi energy. 
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Figure 6: Calculated Fermi surfaces. Top panels: Fermi surfaces across the transition between the weak-coupling SDW and 
the paramagnetic state at conduction band Filing (n) ~ 0.8 and (from left to right) J/t = (0.6; 0.8; 1.0; 1.2; 1.4). The transition 
takes place at J = t. Middle panels: Fermi surfaces inside the SDW phase for (from left to right) Jjt — (0.8; 1.0; 1.2; 1.4) and 
(n) ~ 0.9. The Fermi surface topology changes at J jt — 1.2 (second from right). Bottom panels: Fermi surfaces for different 
dopings across the phase transition between the strong-coupling SDW and the paramagnetic state at J jt — 1.8 for (from left 
to right) (n) = (0.98; 0.93; 0.85). The magnetic order vanishes in the middle panel. 
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Figure 7: First order phase transition between the SDW phase 
and the paramagnetic phase at weak coupling and (n) 0.8. 

The spin polarization as well as the charge density order van¬ 
ish abruptly at Jjt ~ 1. The figure shows the electron po¬ 
larization (top), the spin-spin-correlation, between 

conduction electrons and localized moments (middle), and the 
double occupancy of the conduction electrons (bottom). 


similar change in the energy-momentum dispersion can 
be observed. At weak coupling, the non-interacting 
bands approaching the Fermi energy at (tt, 0), (0, tt), and 
(7r/2,7r/2) are gapped out at the Fermi energy. Bands 
at (tt, tt) and (0,0) do not exist or are very weak. At 
approximately J jt ~ 1.4, this dispersion changes. The 
bands approaching the Fermi energy are bent towards 
each other. The spectral weight at (tt, 0) and (0, tt) close 
to the Fermi energy is shifted to (7r,7r). Although none 
of these bands crosses the Fermi energy, we nevertheless 
observe a similar change in the energy-momentum dis¬ 
persion at half filling as for the doped system. 


VI. PHASE TRANSITIONS 

As we have shown above, the Fermi surface changes 
within the SDW phase. We will show now that the phase 
transition between the SDW phase and the paramagnetic 
phase is also influenced by this change inside the SDW 
phase. In Fig. we analyze the transition between the 
SDW state and the paramagnetic state at weak coupling 
and conduction band filling (n) ~ 0.8. This represents 
an example for the transition between the small Fermi- 
surface SDW and the large Fermi-surface paramagnetic 
state. We show in Fig. the polarization of the con¬ 
duction electrons (top panel), the spin-spin correlation 
between localized spins and conduction electrons (mid¬ 
dle panel), and the double occupancy of the conduction 
electrons (bottom panel). We only show the expecta¬ 
tion values corresponding to the maximally polarized re¬ 
gions of the SDW. As can be clearly seen, the electron 



<n> 

Figure 8: Second order phase transition at J jt — 1.8 for 
different conduction band fillings. The figure shows the same 
quantities as in Fig. 

polarization vanishes abruptly around J jt = 1, which 
signifies a first order transition between the SDW phase 
and the paramagnetic phase. We note that within the 
SDW phase also the electrons in the low electron-density 
regions are polarized within this parameter region and 
that also this polarization jumps abruptly to zero at the 
phase transition. Exactly at this interaction strength, 
the Fermi surface changes from small to large. We also 
observe small discontinuities in the spin-spin correlation 
and the double occupancy at the phase transition. How¬ 
ever, both discontinuities are not very significant. We 
furthermore note that although the Fermi surface be¬ 
comes small, the spin-spin correlation does not vanish 
at weak coupling. 

The transition at strong coupling between the SDW 
phase and the paramagnetic phase looks completely dif¬ 
ferent. We show the polarization of the conduction elec¬ 
trons (top panel), the spin-spin correlation between local¬ 
ized spins and conduction electrons (middle panel), and 
the double occupancy of the conduction electrons (bot¬ 
tom panel) for conduction-band fillings across the phase 
transition at J/t = 1.8 in Fig. Contrary to the phase 
transition at weak coupling, the electron polarization be¬ 
haves smoothly upon doping at strong coupling. We ob¬ 
serve the well-known second order phase transition be¬ 
tween the SDW phase and the paramagnetic phase. Also 
the spin-spin correlation and the double occupancy seem 
to be smooth across the phase transition, although there 
might be small kinks. 

Finally, let us take a look at expectation values across 
the transition of the Fermi surface inside the SDW phase. 
Figure shows the same quantities as shown in Fig. 
and Fig. but for (n) 0.9 and different interaction 

strengths. We have marked the coupling strengths, at 
which the Fermi surface changes from small to large. At 
these interaction strengths, we do not observe any signs 









Figure 9: The same quantities as in Fig. and Fig. 
across the Lifshitz transition within the SDW phase. The 
shaded area corresponds to the interaction strength at which 
the Fermi surface changes from small to large. 


of a transition within the shown expectation values. All 
quantities including polarization as well as spin-spin cor¬ 
relations behave smoothly without any kink. Only for 
large interaction strength, J/t ~ 1.6, we observe the sec¬ 
ond order phase transition between the SDW phase with 
large Fermi surface and the paramagnetic phase. Because 
these results are for the doped system, this transition oc¬ 
curs already for Jjt « 1.6 compared to Jjt ~ 2.2 at half 
filling. 

Summarizing our results for the static expectation val¬ 
ues, we find a first order phase transition between the 
SDW phase with small Fermi surface at weak coupling 
and the paramagnetic phase, and a second order phase 
transition between the SDW phase with large Fermi sur¬ 
face at strong coupling and the paramagnetic phase. The 
change between these two transitions seems to coincide 
with the change in the Fermi-surface topology occurring 
inside the SDW phase. However, we do not observe any 
sign in the static expectation values, although the Fermi 
surface changes inside the SDW phase. 


VII. CONCLUSIONS 

We have analyzed the phase diagram of the Kondo 
lattice model close to half filling by RDMFT. We have 
demonstrated the existence of incommensurate SDW 
states which coexist with CDWs in the Kondo lattice 
model for a wide range of parameters. Although the ex¬ 
istence of SDW states can be expected away from half 
filling, most of the previous calculations have neglected 
this possibility. 

Remarkably, we have found that within this metal¬ 
lic SDW phase the Fermi surface changes from small at 
weak coupling to large at strong coupling. This Lifshitz 


transition, which does not seem to affect static quan¬ 
tities such as polarization and spin-spin correlations, is 
accompanied by a change in the order of the phase tran¬ 
sition between the SDW phase and the heavy fermion 
paramagnetic phase. While we have found a first or¬ 
der transition between the weak-coupling small Fermi- 
surface SDW state and the large Fermi-surface param¬ 
agnetic state, we have found a continuous transition be¬ 
tween the large Fermi-surface SDW at strong coupling 
and the large-Fermi-surface paramagnetic state. 

We want to note that a similar transition between a 
large Fermi-surface and a small Fermi-surface antiferro¬ 
magnetic state has been found in variational Monte Carlo 
(VMC) calculations,I^^H^ although in these calculations 
incommensurate SDW states have not been analyzed. 
These VMC calculations have also shown that the order 
of the phase transition between the paramagnetic state 
and the antiferromagnetic state changes from first order 
at weak coupling to second order at strong coupling. Fur¬ 
thermore, these VMC calculations have confirmed the 
Lifshitz transition inside the antiferromagnetic phase. 
Thus, our RDMFT calculations, which take incommensu¬ 
rate SDW states into account, agree with these previous 
VMC calculations. Antiferromagnetic states with small 
and large Fermi surface have also been found in a differ¬ 
ent DMFT studj^and dynamical cluster approximation 
(DCA) studie#^^^ using the Neel state away from half 
filling instead of SDWs. Finally, we want to note that a 
similar Lifshitz transition has been observed in the fer¬ 
romagnetic phase of the Kondo lattice model.l^ Thus, 
a change in the Fermi-surface topology when increasing 
the interaction strength seems to be a general property 
of the Kondo lattice model and not specific to the SDW 
phase close to half filling. 

Comparing our results to the global phase diagram for 
heavy fermions, we observe two important differences in 
our results. Although we find an SDW phase with small 
Fermi surface and the other with large Fermi surface, 
the phase transition between the SDW phase with small 
Fermi surface and the paramagnetic state is of first or¬ 
der, while it is of second order in the global phase dia¬ 
gram. Second, the spin-spin correlations do not vanish in 
the SDW phase with small Fermi surface. Thus, we do 
not observe a local quantum critical point in the Kondo 
lattice model within the DMFT approximation, which 
is supposed to be between the small Fermi-surface SDW 
phase and the paramagnetic phase. This phase transition 
turns out to be first order in our calculations. We only 
find the previously known continuous quantum phase 
transition at strong coupling between the large Fermi- 
surface SDW and the paramagnetic phase. However, we 
cannot rule out that a local quantum critical point might 
occur in the frustrated KLM, or including long range- 
interactions, or by taking into account spatial fluctua¬ 
tions, which are not included in this DMFT study. By 
including for example an inter-site exchange, the strength 
of the RKKY interaction will be increased, while the 
Kondo effect is not strongly influenced. One can thus 
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expect that the antiferromagnetic phase, and especially 
the small Fermi-surface antiferromagnetic phase, extends 
to larger interactions strengths. Thus, it might be possi¬ 
ble in that case that the Lifshitz transition merges with 
the continuous quantum phase transition. 
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